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We study new physics effects on B decay processes including a final r particle, 
namely B — > Dru and B — > tv. An important feature of these processes is that a 
charged Higgs boson can contribute to the decay amplitude at the tree level in models 
such as Two Higgs Doublet Model and the Minimal Supersymmetric Standard Model 
(MSSM). We derive a resummed effective Lagrangian for charged-Higgs mediated in- 
teractions in the MSSM with the Minimal Flavor Violation. Including supersymmetric 
(SUSY) loop corrections for down-type-quark and charged-lepton Yukawa couplings, 
we calculate the branching ratios of the B — > Dru and B — > tv processes. We find 
that SUSY correction due to gluino-sbottom diagrams can change the Higgs exchange 
contribution by ±50%, whereas stau-neutralino diagrams can make corrections up to 
20%. We also discuss relationship between SUSY corrections in the tauonic decays and 
flavor changing neutral current processes such as B s — > ^ and b — > 57. 
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§1. Introduction 



Recent success of B factory experiments at KEK and SLAC has proved that B physics 
provide a very promising way to explore physics in and beyond the Standard Model (SM). 
The Kobayashi-Maskawa mechanism 1 ) of the CP violation in the quark sector has been 
established from the precise determination of the CP asymmetry in B — > J / ipK$ and related 
modes. 2 )' 3 ) B factory experiments have made many new observations such as the branching 
ratio of the b — > sll 4 ^ and CP violation in the B — > (j)K s mode, 5 )' 6 ) which are known to be 
sensitive to new physics effects. 

In future, more information on B decays will be obtained at current B factories as well as 
hadron B experiments at Tevatron and LHC. Furthermore, the future upgrade of the e + e~ 
asymmetric B factory, Super B Factory, is discussed, where the goal of the luminosity is 
50-100 times more than the current achieved luminosity. 7 ) 

We study here new physics effects on B decay processes including a final r particle, 
namely B — > Dtv and B — > tv. An important feature of these processes is that a charged 
Higgs boson can contribute to the decay amplitude at the tree level in models such as Two 
Higgs Doublet Model (2HDM) and the Minimal Supersymmetric Standard Model (MSSM). 
From the experimental side, since at least two neutrinos are present in the final state in 
the signal side, full-reconstruction is required for the B decay in the opposite side. For the 
B — > Dtv process, the branching ratio is not yet measured even though the SM prediction 
is 8 x 10 -3 . (The inclusive b — > crv branching ratio was determined at LEP experiments. 8 )) 
The B — > tv process has a smaller branching ratio (~ 9 x 10~ 5 ) in the SM, because of 
the helicity suppression, and the upper bound of 2.9 x 10~ 4 is reported from the BELLE 
experiment 9 ) and 3.3 x 10~ 4 from the BaBar. 10 ) These processes will be important targets 
of coming B factory experiments. 

In this paper, we calculate the branching ratio of B — > Dtv and B — > tv processes 
in the MSSM, taking account of supersymmetric (SUSY) corrections to the charged Higgs 
Yukawa couplings. At the tree level, the Higgs sector of the MSSM is of the same form as 
the type II 2HDM, where one Higgs doublet provides mass terms for up-type quarks, and 
the other does for down-type quarks and charged leptons. SUSY loop corrections, however, 
can induce Yukawa couplings of the opposite type. 11 ) -13 ) In particular, it is known that 
new contributions to the bottom and tau Yukawa coupling constants induce flavor changing 
processes such as B s — > /i + /i~, 14 )' 15 ) b — > s/Z, 16 ) r — > 3/Lt, 17 )' 18 ) r — > fir/ 19 ^ and the fx-e 
conversion at muonic atoms, 20 ) especially for a large ratio of two vacuum expectation values 
(tan/?). For the b — > c(u) transition, SUSY corrections were investigated in the inclusive 
b — > ctz/ 21 ) and the B — > tv processes. 22 ) In this paper, in addition to the correction to 
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the charged Higgs boson(iJ ± )-6-c('u) vertex, we include the correction to the H ± -t-v vertex, 
and study importance of these corrections to these two processes. As for the flavor mixing 
in the squark and slepton sectors, we take the assumption of the Minimal Flavor Violation 
(MFV), where a unique origin of the flavor mixings is given by the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix *) . We show that the corrections to the B — > Drv and B — > rv 
branching ratios are large for tan/5 > 30. The tau vertex correction can give sizable effects 
for reasonable parameter sets of squark and sleptons. We also consider correlation of these 
processes to b — > 57 and B s — > fi + fi~ . 

This paper is organized as follows. In section 2, we derive effective Yukawa interactions 
of the charged Higgs boson and quarks/leptons taking account of SUSY loop corrections in 
the MSSM. We apply this formulation to calculate the B — > Drv and B — > tv branching 
ratios. In section 3, results of numerical calculations are presented including the b — > S7 and 
B s — > fi + fi~ . Conclusions are given in section 4. Derivation of the charged Higgs coupling 
in the MFV case is given in Appendix A. 

§2. Formalism 

2.1. SUSY loop corrections to Yukawa interactions 

In this section, we derive the general form of the effective Lagrangian for the charged 
Higgs Yukawa coupling in the MSSM taking account of SUSY loop diagrams. For the flavor 
mixing in the squark sector, we take a model based on the assumption of MFV, where the 
CKM matrix is the only source of flavor and CP violations. As we see below, there is a new 
contribution to the Yukawa coupling constant in a large tan f3 regime. We resum the tan (3- 
enhanced contributions consistently for the down-type and charged-lepton Yukawa coupling 
constants following the method developed in the reference by A.Dedes and A.Pilaftsis. 15 ) 

The superpotential of the model is given by 



W = -H.D^Q + H 2 U c y u Q - H 1 E c y e L + nH x H 2 , (2-1) 
where the components of weak doublet fields are denoted as 




The quantum numbers of SU(3) x SU(2) x U(l) gauge groups for H ± , H 2 , Q, L, D c , U c , E c 
are (1,2,-1), (1,2,1), (3, 2,|), (1,2,-1), (3, 1, §), (3,1,-f), (1, 1, 2). The soft SUSY 
*) For general SUSY models, the tauonic B decay is considered in G.H.Wu et al. 23 ^ 
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breaking masses and trilinear SUSY breaking terms (A-term) are given by 

Aoft = -Q[M| L g L - UlMlJj R - DlMl R D R - Z[m| l Z l - E R M 2 ~J R 

+HiD R A d Q L - H 2 U R A U Q L + H^Aja, + h.c, (2-3) 

where the fields with tilde (~) denote squarks and sleptons. 

In this section, we first discuss the simplest case where soft breaking masses are propor- 
tional to a unit matrix in the flavor space, and A u , A d and A e are proportional to Yukawa 
couplings. Their explicit forms are given below, 
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(2-9) 
(2-10) 
(2-11) 

where a,i(i — 1 — 5) are real parameters. 

At the tree level, the Yukawa couplings have the same structure as the above superpo- 
tential, namely, Hi couples to D c and E c , and H 2 to U c . On the other hand, different types 
of couplings are induced when we take into account SUSY breaking effects through one-loop 
diagrams. Lagrangian of the Yukawa sector can be written as 

^Yukawa = —H 1 D R y d Q L + H 2 U R y u Q L — HiE R y e L L 

-H 2 D R Ay d Q L + HiU R Ay u Q L - H 2 E R Ay e L L + h.c. (2-12) 
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where Hip = icr 2 Hl 2 , and Ay d , Ay u , and Ay e are one-loop induced coupling constants. 
(Here and in the followings, quark and lepton fields in capital letters represent three vectors 
in the flavor space.) 

From the above Yukawa couplings, we can derive the quark and lepton mass matrices 
and their charged Higgs couplings. For the quark sector, we get 

£quark = cos (5D R y d [l + tm/3A md ]D L + smf3H~D R y d [l - cot f3A md ]U L 

--j= sin/3[7 R y u [l - cot (3A mu ]U L + cos/?# + [7 R y u [l + tan/3Anj£> L + h - c -> 

(2-13) 

where we define A md (A mu ) as A„ ld = y^Aya (A rriu = y^Ay u ), and v ~ 246GeV. Notice 
that Ay d is proportional to or y^y^y^ in this case. We then rotate the quark bases as 
follows: 

U L = V L (Q)U L ', D L = V L (Q)V CKM D L ', (2-14) 
U R = V R (U)U R , D R = V R (D)D' R , (2-15) 

where the fields with a prime ( ' ) are mass eigenstates. In this basis, the down-type quark 
Lagrangian is given by 

AD-quark = cos PD R vl(D)y d V L (Q)R d V C KMD L ' 

+ smf3H~D R Vl(D)y d V L (Q)U L ' + h.c, (2-16) 

where R d = 1 + tan/3 Z\ md and A mi = V^(Q)A md Vi J (Q). Hereafter, a matrix with a hat (") 
represents a diagonal matrix. Since the down-type diagonal mass term is given by 

M d = ^=cos^(D)y d V h (Q)R d VcKM, (2-17) 



we obtain the following Lagrangian for down-type quarks. 

— /2 

£ D -quark = -D R M d D L ' + — tan pH-D^MaV^RfU^ + h.c. (2-18) 



v 

The corresponding corrections to the up- type couplings can be calculated from Eq.(2-13). 
Since we are interested in the large tan (3 case, these corrections are very small. In the 
following, we neglect such corrections, and the Lagrangian for the up-type-quarks is given 
as follows: 

£u-quark = -U R ' M V U L ' + ^- cot (3 H + U R M U V CKM D L ' + h.c. (2-19) 
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For the case of the charged-lepton, we can derive relevant parts of the Lagrangian in a 
similar way to the case of the down-type quark. Using the following definitions, 



E L = V L (L)E' L , E R = V R (E)E' R , (2-20) 
A me = y~ e l Ay e , A nie = V^(L)A me V L (L), (2-21) 

we obtain the Lagrangian for the charged lepton as follows: 

/2 

Acpton = -E' R M e E' L + y- tan (3H~ e' r M e R~ l N L + h.c. . (2-22) 

Here, the prime represents the mass eigenstate, and we neglect the neutrino masses, and 
R e = 1 +tan/3Z\ me . 

In the present case with Eq.(2-4) - Eq.(2-10), A md receives contributions from gluino and 
down-type squark, and higgsino and up-type squark diagrams. The explicit form is given as 
follows: 



E~ g + E~ h , (2-23) 



where 



2a 

E~ = ^l^M g /[Mg, M 3l , M 5 J, (2-24) 
\ = -T^r 2 A ^\ 2 I[M~ h ,M Dij ,M D J, (2-25) 



a 2 6 2 lng + 6 2 c 2 ln^ + c 2 a 2 ln4 , . 

^ - ^ W -^-^) a - (2 ' 26) 

and E^ are gluino and charged higgsino contributions shown in Fig. 1(a) and (b), re- 
spectively. Note that these corrections for Yukawa couplings are calculated in the unbroken 
phase of SU(2) x U(l) symmetry. For the charged lepton case, A m& is given by 



An e = E 5 = {M \ v2 ^ ] ^M~ B I[M~ B , M £l , M £ J, (2-27) 

from the bino-slepton diagram shown in Fig. 2. 

Up to now, we have assumed all squark mass matrices are proportional to a unit matrix 
at the electro- weak scale, as shown in Eq.(2-4) - Eq.(2-8). However models with the MFV 
correspond to more general cases. For instance, the assumption of Eq.(2-4) - Eq.(2-8) is not 
satisfied in the minimal supergravity where all squarks have a universal mass at the Planck 



6 



Dr. 



9r 9l 



H° 2 * 



(a) (b) 

Fig. 1. Non-holomorphic radiative corrections to the down- type quark Yukawa couplings induced 

by (a) gluino <?l,r and (b) charged higgsino h\ 2 . 
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Fig. 2. Non-holomorphic radiative corrections to the charged lepton Yukawa couplings induced by 
bino -E>l,r- 

scale, not at the electro-weak scale. In appendix we derive the charged Higgs coupling in 
more general case of MFV. Namely the squark mass matrix is taken to be 

M| l = [ail + b iy ly u + b 2 y\y d ]M\ (2-28) 
M 2 ~ n = [a 2 l + hy u y{]M\ (2-29) 
I 2 ~ 

D R 

The final results of the charged Higgs coupling is given by 



m = [a 3 l + b 6 y d y d ]M 2 . (2-30) 



C H± *V-ttmf3H-D R ' ^ * Vi KMij U L ' j + h.c. for (ij) = (1, 1), (1, 2), (2, 1), (2, 2), 



V2 + M di 
— tan pH D Ri p- 

v H *1 + [£ 5 w ]tan/3 

(2-31) 

C H± « ^ tan/?g-D R ; i - [e/) Md EL{i])] — 4m,^ + h.c. for (i,,) = (3,1), (3, 2), 

(2-32) 

r )H + RM-n 1 i + [ jg J3) + ^(33 )]tan/ g 

* t; ^ Rl 1 + tan/? 1 + [EM + + ^ W) + ^.(«) + (as)] tan/ j V ckm^l 3 

+h.c. for (i,;) = (1,3), (2, 3), (2-33) 
£„ ± « ^ tan/Mf-Z^ — ^ ^^-— ^^ + h.c. for = (3, 3). (2-34) 
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where the function Eqf l \ etc are listed in Appendix A. In deriving these results we only 
keep yt in the up-type Yukawa coupling in loop diagrams and use the hierarchy of the CKM 
matrix elements. See Appendix A for details. Notice that the above results do not depend 
on the relation between the A-terms and the Yukawa couplings, since we only keep y t in loop 
diagrams, even though Eqs.(2-9) and (2-10) are assumed in Appendix A. 

2.2. Effective Lagrangian for J^-Fermi interactions 

Once we obtain the charged Higgs and fermion couplings, it is straightforward to write 
down the amplitudes for B — > Dtv (B~ — > D t~v or B — > D + t~v) and B — > tv processes. 
First, the effective Lagrangian for b — > ctv operators is given by 

G F _ 

£ c s = -—^VcbC-f^l - 75)6x7^(1 - 7 5 )z/ T 

+G s cbT(l - 75 )z/ r + GpC7 5 &r(l - l5 )v T + h.c, (2-35) 

where G$ and Gp are scalar and pseudo-scalar effective couplings. These couplings are given 
from Eqs.(2-18), (2-19) and (2-22), 

G s = t -^^[K%(M b [R- d 1 h2V cb + M c V cfe cot 2 /3), (2-36) 

G P = A%3(M b [R^] 22 V cb - M c V cb cot 2 (3). (2-37) 

We omit a prime ( ' ) from the fields in mass eigenstates. The higgsino diagram contributions 
to the [-R^" 1 ]22 is proportional to square of the charm Yukawa couplings. Since the branching 
ratio can change only by at most a few %, we neglect such contributions in the followings. 
For the case of large tan/3, we can also neglect the last terms in Gs and Gp. 

In order to calculate the B — > Dtv branching ratio, we need vector and scalar form factors 
of the B — > D transition. In the heavy quark limit, these form factors can be parameterized 
by a unique function called the Isgur-Wise function. The form of the Isgur-Wise function was 
investigated by using the dispersion relation. 24 ) From the semi-leptonic decays B — > Dlv 
and B — > D*lv (I = e, //), the Isgur-Wise function is obtained in a one-parameter form, 
including the short distance and 1/Mq (Q = b, c) corrections. The short distance corrections 
for B — > Dtv was also calculated. 25 ^ Here we adopt the Isgur-Wise function obtained in 
these literatures, but we do not include the short distance and the 1/Mq corrections for 
simplicity. The short distance effects were shown to change the branching ratio within 6% 
in the reference of T.Miki et al. 25 ) 

Using the definitions, 

2p B .p D _ 2 PB .p T Ml M 2 
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the differential decay width is given by 



^ = -128^ Mb ^^^' (2 ' 39) 

where 

p D (x,y) = [\f + \ 2 9l (x,y) + 2Re(f + f*)g 2 (x,y) + \fL\ 2 g 3 (x)i (240) 
gi(x,y) = (3-x-2y- r D + r T )(x + 2y - 1 - r D - r T ) 

-(1 +x + r D )(l + r D -r r -x), (2-41) 

g 2 (x,y) = r T (3 - x - 2y - r D + r T ), (2-42) 

030*0 = r T (l + r D - r r - x), (2-43) 

/: = [/_ - A S [U(1 - r D ) +/_(! + rD - *)]], (2-44) 

/± = ± i±^ e(w ) ; (2-45) 

Here As = GfV bM ( M ^_ M ) • We use the following form of the Isgur-Wise function. 

£(w) = 1 - % p \z + (51/o? - 10)z 2 - (252p 2 _ 84)^ 3 , (2-46) 

t = CT-^ (2 . 47) 

For the slope parameter, we use p 2 = 1.33 ± 0.22. 24 )' 25 ) 

For the B — > rz/ process, the relevant four fermion interactions are those of the 6 — > -urz/ 

type, 



G 

£' cS = -^|KbM7 M (l - 75)6x7^(1 - 7 5 )z/ T 

+G" s U6r(l - 7 5 )z/ r + GpU756r(l - 7 5 )z/ r + h.c, (2-48) 
G' s = T fe 1 ] 33 (M fe fe 1 ] 11 K fc + M u K 6 cot 2 /5), (2-49) 

G"p = [4" 1 ]33(M 6 [^- 1 ] 11 K 6 - M u V ub cot 2 {3). (2-50) 



Using the matrix elements 



(0\u^ 5 b\B-) = if B p», (2-51) 

M 2 

(0\u l5 b\B~) = -ifBjf, (2-52) 
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the decay width is given by 



G 



2 



r[B - rv\ = ^\V ub \ 2 flMlM 2 T 



2 



i 2v n> 

M b M T V ub p 



2 



1-rv) 2 , (2-53) 



where /b is the _B U decay constant. 

In the generalized case of the MFV with Eqs(2-28) - (2-30), the scalar and pseudo-scalar 
couplings, Eqs.(2-36), (2-37), (2-49), and (2-50) can be obtained by the following replacement. 

[R ~ dl]22 ~* l + [E/V -ELW]t^ (2 ' 54) 

[R " 1]n ~* 1 + [Es® -ELW] tan P (2 ' 55) 

Notice that the right-handed sides of the above equations are approximately same because 
fts E~- 32 \ This is the generalization of [-R^^n ~ [-R^ 1 ^, which follows from fact that 
the higgsino diagram contribution can be neglected in the evaluation with the [-R^n and 

§3. Numerical results 

In this section, we shall present results of the numerical calculations on branching ratios 
of the B — > Dtv and B — > tv processes in the MSSM. We see that charged Higgs effects 
to these processes become important for the parameter region of a large tan (3 and a small 
charged Higgs mass. We also discuss B s — ■> and b — > 57, because SUSY corrections to 

these processes are important for this parameter region. The relevant SUSY parameters are 
tan/5, Mh±, the higgsino mass parameter /x, the bino mass parameter Mg, the gluino mass 
Mg, the sbottom mass Mr, and the stau mass M?. For the sbottom and the stau, we take 
the left and right handed masses to be the same, and neglect the left-right mixing terms. 

We first show [R^ 1 ^ i n Eq.(2-18) and [i?^ 1 ^ in Eq.(2-22). As we discussed in the 
previous section, [-R^" 1 ]].! and [-R^ 1 ^ are approximately same, because the higgsino loop 
contributions, Eq.(2-25), are suppressed. Therefore the A-term dependence of [R^ 1 ] 11,22 is 
negligible. In Fig 3, contour plots of [-R^ 1 ^ are presented in the gluino mass and light 
sbottom mass eigenvalue space in the case of the \i = ±400GeV and tan/5 = 50. For a 
positive value of /i, the correction become —10 to —40 % in this parameter region. On the 
other hand, corrections become positive and huge for a negative ji. The value of the [-R^T 1 ^ 
is shown in the light stau mass eigenvalue M ?1 and bino mass space for /x = ±400GeV and 
tan/3 = 50 in Fig. 4. In general, the correction is smaller compared with the case of [-R^ 1 ^- 
For a larger value of /x, however the correction can be 10 %. One such example is shown in 
Fig.5. 
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fj, = 400GeV, tan/5 = 50 
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Fig. 3. Contour plots of the correction factor [R^^uii = 1, 2) in the Mg and space for tan f3 = 
50 and fj, = ±400GeV. The numbers in figure are the values of [K^ l }n(i = 1, 2). The values of 
[R^^uii = 1,2) is unity without SUSY corrections. 

The dependence on tan/3 in and [-R" 1 ^ are shown Fig.6. Here we take /x = 

±400GeV and /i = ±200GeV. The gluino and bino masses satisfy the GUT relation (Mg = 
Q.72Mg). The correction to [R d 1 }22 becomes large for tan (3 > 30. 

Next, we show the branching ratio of B — > Dtu and B — > rv. For B — > Drz/, we 
consider the following ratio of the two branching ratios, where B(B — > D/ji/) is denned as 
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Fig. 4. Contour plots of the correction factor [-R^T 1 ^ in the and for tan/3 = 50 and 
fi = ±400GeV. The numbers in figure are the values of [i?" 1 ^. The value of [i?^ 1 ]33 is unity 
without SUSY corrections. 



the branching fraction for the B — > Dfiu mode integrated over the same phase space of the 
B — > Dtv kinematics. 



B(B -> Dry) _ f]%T dy^M^y) 
B(B - Dfiv) dx ff dy$^M*f/ D (x, y) ' 
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Fig. 5. Contour plots as Fig.4 for tan/3 = 50 and \i 
the bino and light stau mass eigenvalue. 



±2000GeV. We also change the range of 



J/1,2 = 



1 - I ± 



(3-2) 
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Fig. 6. The correction factors as a function of tan/3. The upper figure is [R d ]u(i = 1,2) for 
Mg = 672GeV and M~ b = 400GeV, and the lower figure is [R-%3 for M § = lOOGeV and 
M^f = 300GeV. The value of the fi parameter is shown in the figures. 
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For the B — > D/iz/ mode, we use 



Pkfoy) = [IZ+IYifol/) + 2Re(f + f*)g' 2 (x,y) + l/ilVaW], 



(3-4) 
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g[(x,y) = (3-x-2y- r D + r^)(x + 2y - 1 - r D - r„) 

-(1 +x + r D )(l +r D - r M - x), (3-5) 

g' 2 ( x ^y) = r fJ ,(S-x-2y-r B + r fJ ), (3-6) 

03 (a) = r^l + rn -r^ - x), (3-7) 



M, 



It was pointed out that we can reduce theoretical uncertainty associated with form factors by 
taking this ratio. 26 ) In Fig. 7, we show the above quantity as a function of the charged Higgs 
mass for the parameter set of M~ h = 400GeV, M § = lOOGeV, M 7 = 300GeV and tan 3 = 30 
and 50. In this figure, we also draw a line without the SUSY corrections. We can see a 
large deviation due to the SUSY loop corrections. A similar figure for B — > tv are shown 
in Fig.8. We can see that B(B — > tv) vanishes in a particular point of M H ± depending on 
SUSY parameters, and below that point the branching ratio is significantly enhanced. 

We consider the correlation between the branching ratio of the B — > Dtv and that of 
the B — > tv. Under the assumption of MFV, the charged Higgs effect appears through the 
following combination of the parameters in the branching ratio formulas, 

^^^V^'bfcV (3-9) 

(or the replacement of [B~ L 1 } 2 2 by the right hand side of Eq.(2-54) in the generalized MFV 
case.) This is compared with the type II 2HDM, where R is replaced by R, 

M w tan 3 . 

R = ^k^- (3 ' 10) 

In other words, the SUSY corrections effectively change the value of the tan 3 in the formula 
of the 2HDM. Therefore the correlation between two branching ratios is the same for MSSM 
and 2HDM. Note that this correlation arises due to our assumption of the MFV, namely 
from the fact that the higgsino contribution does not induce sizable effects in these branching 
ratios. In Fig. 9, we show the correlation of the two quantities. We also show estimated 
theoretical uncertainties for several values of R along the line. The error corresponds to 
uncertainty from p\ = 1.33 ± 0.22 for B -> Dtv, and from f B = 200 ± 30MeV and \V ub \ = 
(3.67 ± 0.47) x 10~ 3,27) for B -> tv. We see that uncertainty of B(B -> tv) from the 
present input parameters is still large. These errors, however, can be reduced significantly 
in future from more precise determination of semi-leptonic B decay form factors and \V u b\, 
and improvement of determination from the lattice gauge theory. For instance, if the 
uncertainty of the slope parameter p\ is improved to ±0.10, the theoretical uncertainty from 
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Fig. 7. Value of ^ b ^ Dtu ) as a function of the charged Higgs mass for tan/? = 30,50, Ms 

B(B •/»;/;■• 



lOOGeV, Mt; 



6.72M 5 , 



400GeV and M= = 300GeV. The horizontal solid line is the 



predicted value in the SM. 



this source to B(B — > Dtu) is reduced from ±5% to ±2% for the SM case. If the uncertainty 
of \V u b\ is improved to ±5.8%, the error of B(B — > tu) is reduced from ±39% to ±32%, and 
further improvement of \V ub \ and f B to ±4.4% and ±16MeV respectively leads to ±18%. ^ 
At the SuperKEKB, it is expected that the sensitivity to R reaches to 11 (90% confidence 
level) for an integrated luminosity of 5ab _1 from the B — > Drv process. 7 ) Observation of 
the B — ► tv mode is possible at 30ab _1 for the SM case. 28 ) 
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M s = lOOGeV, Mg = 6.72M 5 , M~ b = 400GeV, Mj = 300GeV and f B = 200MeV. The 
horizontal solid line is the predicted value in the SM. 



Let us comment on scaling behaviors of the SUSY loop corrections. E^, E^, and E^ in 
Eqs.(2-24), (2-25), and (2-27) remain constants when all SUSY mass parameters are multi- 
plied by a same factor. Therefore the SUSY loop effects to the charged Higgs contribution in 
B — > Dtu and B — > tv do not decouple by taking large SUSY mass spectrum as long as the 
charged Higgs mass is the same. This situation is similar to the SUSY loop contributions to 
the neutral Higgs exchange in the flavor changing neutral current and lepton flavor violation 
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B — > Dtu/B — ► .D//^ and and |V^&| for B ^ tv are also shown for each point. 



processes. 14 ^ 20 ) 

Finally, we discuss correlation of the tauonic B decays with b — > S7 and £> s — > fi + fi~. 
It is known that these processes receive significant SUSY contributions for the large tan/5 
case. In the MFV, the b — > 37 amplitude consists of the SM contribution, the charged Higgs 
contribution, and the chargino-stop contribution. The gluino-sbottom contribution is not 
significant for MFV. The effect of SUSY loop correction to the charged Higgs vertex was also 
studied. 29 )' 30 ) We calculate the b — * S7 branching ratios following the formula presented in 
G.Degrassi et al. 29 ) Since the b — * 57 process depends on the chargino-stop diagram, there is 
no strict correlation between b — > 57 and tauonic B decays. However, the contribution from 
the charged Higgs diagrams is enhanced for fi < from the correction by R~ l in Eq.(2-18). 
As an example, we show B(b — > 57) for the following parameter sets, = lOOGeV, the wino 



mass Mr, 



1.99Mg, Mg 



6.72M 5 , M~ t 



M~ b = 400GeV, \A U \ = lOOGeV, |/x| = 400GeV, 



and tan/5 = 50 in Fig. 10. In order to satisfy the experimental constraint, a fine tuning 
between the charged Higgs and chargino contributions is necessary for /j, < 0. For /i > 0, 
the constraint is generally weak. 

We also calculate the branching ratio of the B s — > in this model. In this case, 

the branching ratio depends on the higgsino-stop diagram and the gluino-sbottom diagram. 
In particular, inclusion of the higgsino loop contribution is necessary to generate this fla- 
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Fig. 10. b — > S7 branching ratio as a function of the charged Higgs mass for = lOOGeV, 
M w = 1.99Mg, Mg = 6.72M § , M~ t = M~ b = 400GeV, \A U \ = lOOGeV, \n\ = 400GeV, and 
tan/3 = 50. The shaded region is the experimental allowed region at 2a level. 31 ) 

vor changing process. We present the branching ratio of B(B S — > in Fig. 11 for the 

parameter sets, tan/3 = 50, \A U \ = \A d \ = lOOGeV, = 400GeV, M~ t = M~ b = 400GeV, 
M § = lOOGeV, Mg = 6.72M 5 , = 300GeV, and f Bs = 230MeV. The present experimen- 
tal upper bound, 7.5 x 10 -7 ', 32 ) is also shown. We can see that the large branching ratio is 
expected for the parameter space where the tauonic B decays receive significant contribution 
from the charged Higgs diagram. The enhancement is particularly large for \x < 0. 
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Fig. 11. Branching ratio of the B s — > as a function of the charged Higgs mass for tan (3 = 50, 

\A U \ = \A d \ = lOOGeV, |/i| = 400GeV, M~ t = M~ b = 400GeV, Mg = lOOGeV, M~ g = 6.72Mg, 
Mji = 300GeV, and fs s = 230MeV. The horizontal solid line is the current experimental upper 
bound. 32 ) 

The SUSY loop corrections can be also important in the anomalous magnetic moment of 
muons. In this case the relevant SUSY parameters are the slepton masses, gaugino masses, //, 
and tan (3. For example, the parameter set taken in Fig. 10 and 11, the SUSY contributions to 
(g — 2)^/2 is about 7 x 10 -9 , which is larger than the current discrepancy between the exper- 
imental result and SM prediction. 33 ) If the slepton mass is taken to be larger than 400GeV, 
the SUSY contributions to (g — 2)^/2 can be compatible with the current discrepancy. Notice 
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that slepton contribution to the B s — > and tauonic B decays are sub-dominant and 

that for b — > S7 is negligible. Therefore the SUSY contributions to 5 S — > and tauonic 

B decays can be important while both (g — 2) M /2 and 6 — > S7 constraints are satisfied. 

§4. Conclusions 

In this paper we study SUSY effects on the tauonic B decays, -B — > -Drz/ and 5 — > rz/ 
under the assumption of the MFV. These processes receive large corrections for large tan (3 
regime through SUSY loop diagrams related to the bottom and tau Yukawa couplings. For 
the bottom Yukawa coupling, we find that only the gluino-sbottom loop can contribute 
significantly and the chargino loop contribution is at a few % level. This is in contrast to 
the case of the correction to the b-t-H vertex, where the stop-chargino loop correction is 
also important. 13 ^ The effect of the stau-neutralino loop to the tau vertex is generally not as 
significant as the sbottom-gluino loop, but in some parameter space these contributions can 
change the charged Higgs exchange effect by more than 10 %. We also study the correlation 
between B(B -> Dtv) and B(B -> ru) within the assumption of the MFV. The SUSY 
effect on these processes can be absorbed as an effective change of the tan (3 value, so that 
the correlation itself is the same as the 2HDM without SUSY loops. It would be therefore 
interesting to compare this value with tan j3 measured from other processes in MSSM such 
as the heavy Higgs direct production, 34 ) the chargino/neutralino mixing, 35 ) and the stau 
decay. 36 ) SUSY corrections to the Higgs couplings can also change the light-Higgs branching 
ratios 12 ^ and the search limit of the SUSY Higgs bosons at LHC. From combined analysis 
of Super B Factory and collider experiments, we may be able to obtain important insight 
on flavor mixing for squarks and sleptons. For example in the present case of MFV, SUSY 
corrections to b-c-H and b-u-H verteces are the same while that to the b-t-H vertex can be 
different because of the higgsino loop contribution. 

We also compare the tauonic B decay fractions with b — > 37 and B s — > branching 
ratios. In general, these processes receive large corrections when we expect large effects 
in tauonic B decays, i.e. a large tan/5 and small M H ± region. Since stop and chargino 
diagrams are essential in these flavor changing neutral current processes, we do not have a 
strict correlation among these processes. However the parameter space with ji < is strongly 
constrained by b — > and B s — > . 

The tauonic B decay processes considered here provide important information on the 
Yukawa interaction associated with the charged Higgs boson. A large deviation from the 
SM prediction is expected for large tan (3 cases from both tree level and SUSY loop effects. 
The tauonic B decays at future B factory experiments therefore can play a unique role in 
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exploring SUSY models. 
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Appendix A 

Charged Higgs coupling for the general MFV case 



In section 2, we have derived the resummed effective Lagrangian for the charged Higgs 
couplings under the assumption that soft SUSY breaking mass matrices for the squarks are 
proportional to a unit matrix in the flavor space. In this appendix, we relax this assump- 
tion and obtain the charged Higgs couplings in a more general case of the MFV. General 
consideration on the MFV case is given in the literature. 37 ^' 38 ^ 

We consider the following squark mass matrices and A-terms. This form is motivated 
from the renormalization group effects on the squark mass matrices in the SUSY breaking 
scenarios such as minimal supergravity, gauge mediation, and anomaly mediation. 

M| l = [ail + &iyty u + b 2 y\y d ]M\ (A-l) 

M2 R = [a 2 l + 6 3 y u yi]M 2 , (A2) 

Ml ii = [a 3 l + hy d y\]M 2 , (A3) 

A-uij A u y u ij, (A-4) 

A*,- = A d y dij , (A- 5) 

where 01,2,3 and 61,2,3,4 are real parameters. In the following we take the basis where the 
down-type Yukawa coupling is diagonal. 

Yu = y«u c ° KM , (A-6) 
Yd = Yd, (A- 7) 

where V^km * s the flavor mixing matrix in the original Yukawa coupling. Furthermore when 
we calculate the loop diagrams, we use approximation that only the top Yukawa coupling is 
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Fig. 12. Subleading corrections to the down-type quark Yukawa couplings induced by (a) gluino 
5l,r and (b) charged higgsino h^ 2 - 



kept in the up-type Yukawa coupling. 



/ 



V 



(A-8) 



Vt I 



Explicit form of the mass matrix is given by 

f [ ai + b 2 yj}M 2 ^ 
[a l + b 2 y 2 ]M 2 



M% 



fc»3 



V 



+ 



(A-9) 



[ai + hy 2 + b 2 y 2 b ]M 2 J 

^ll/i^CKMl3^CKM32^ 2 ^ 2 ^CKM13^CKM33^ 2 \ 



°lVt l/ CKM23 K CKM31^ °lVt "CKM23 V CKM 33 IV1 

\ °Wt K CKM33 K CKM31 iW °Wt ^CKMSS K CKM32 iW 



/ 



M .. + AM 2 ~ , 



M 2 ~ 



( a 2 M 2 



\ 



a 2 M 2 



= M 2 ~ . 

Ub.i 



(A-10) 
(A-ll) 



M| 



(A-12) 



[a 2 + 6 3 % 2 ]M 2 / 

/ [a 3 + 6 4 ^]M 2 \ 

[a 3 + b,y 2 ]M 2 
V [a 3 + b A y 2 }M 2 J 

We calculate the correction to the down-type Yukawa coupling in Eq.(2-12). In addition 
to the Fig.l, there are extra contributions shown in Fig. 12, where the off-diagonal terms in 
Eq.(A-lO) are treated as mass insertion. The explicit form of Ay d is given by 

/ Ay du Ay dl2 Ay dl3 \ 



Ay dij 



Ay d21 Ay d22 Ay d23 
\ Ay d31 Ay d32 Ay d33 ) 



(A-13) 



23 



where, 



Ay <m 


= y d EgW 


(A-14) 


Ay dl2 


= y d [EL {12) + E X ( 13 ) + 4 (133) ]^ M13 \/ C ° KM32 , 


(A-15) 




= y d [ELM + + ^ (133) ]\/ C % 13 \/ C KM 33, 


(A-16) 


Ay d 2i 


ee y s [4 (21) + E~W + 4 (233) ]^ M23 \/ C KM31 , 


(A-17) 


Ay d22 




(A-18) 


Ayd23 


— J/sPj + ^ft ^ ft J V CKM23 V CKM33' 


(A-19) 


Ay d3 i 


- „ TP /(31) _L /? ( 33 )l T/°t T/° 

= yb[^g -r & h J ^CKM33 K CKM31> 


(A-20) 


Ay d3 2 


_ „, r 77-/ (32) p (33)i T/ 0t T/° 

= yb[Eg + E^ J* / CKM33 l/ CKM32' 


(A-21) 


Ay d33 


= y b [E~ g (3) + ^ (33) ]^ M 33^C°KM33, 


(A-22) 



h 



2a s n* 
3tt Ms 



j(3) 



M n i M n i 



V A u 2r(3) 
167T 2 Mr y * 



Mr, 



Mr ' Mr 



/ (u ) 



lots \i* 
' 3tt M| 



b lV 2 M 2 I^ 



Ql« Qu £> r « 



E 



i {U3) 



V A u 

16tt 2 M 3 

h 



M~ g Mg 



M Un3 



I<®[x,y] 
lW[x,y,z] 



x 2 In x 2 



y 2 In y 2 



(x 2 



l)(x 2 -y 2 ) 
x 2 In x 2 



(y 2 



i)( y 2 -x 2 ) 



(x 2 — l)(:r 2 — y 2 )(x 2 — z 2 ) 
y 2 \ny 2 



+ 



In z 2 



(A-23) 

(A-24) 

(A-25) 

(A-26) 
(A-27) 



(y 2 - l)(y 2 - x 2 )(y 2 - z 2 ) (z 2 - l)(z 2 - x 2 )(z 2 - y 2 ) 
From the Ay d , we obtain the mass term for up-type and down-type quarks. 



£ qU ark = -D Ri [M dij + A i:j ]D hj - U Ri M uij U Lj + h.c, 



(A-28) 



(A-29) 



here 



M di , = 



( Mn 



V 



M- 
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M- 



33 



/ 



M di , 



M dll = cos/3y d [l + tan/5], 
V2 



(A-30) 
(A-31) 
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M d22 = -j= cos/3y s [l + E~ g M tan/3], 



M d33 = 



cos/3^1 + (£/ 3) +^ (33) )tan/3], 
A, \ 



A. 



v 

f A12 AJ 13 

Al A23 

\ Ai A 32 

(12) 



/ 



cos (3y d [E^+E~ h ^ 



A 12 

A 13 = JL cos «4 (13) +^ 13 
Ai 

A 3 

An 

A 2 



v/2 

f 

71' 

-^cos^[< 1 ) + ^ 23 

-^cos/3y s [4 (23) +^ (23 
-^cos^[4 (31) +^ (33 
-^cos^[4 (32) +^ (33 

71' 



+ -Er^ ^CKM13^CKM32 tan/3, 



, (133) lT/ 0f 



^ckmi3^ckm33 tan/3, 



/ (233) lT/ 0f 



ck \J23^CKM3i tan/3, 



w(233) 1I/0 t t^O f o 

+ jC/ fc J l/ CKM23 K CKM33 Xan / :> ' 

]^ / CKM33^ / CKM3i tan/3, 

]^ / CKM33^ / CKM32 tan/3, 



M u = — sin/3y u V^ KM . 



Next, we rotate the basis to mass eigen-states as follows: 

u L = v L (u)u L ' = v£l M u h ', U R = V R (U)U R ' = U R ', 

D L = V L (D)D L ', D R = V R (D)D' R . 
We introduce AVl and AV R as 

V L (D) = 1 + AV L , 
V R (D) = 1 + AV R , 

where the unitarity requires that ZWjj^ = — AVl(r) and AVl( R )u = 0. At the 
Aij, AVl, r is expressed as 39 ) 



(A-32) 
(A-33) 

(A-34) 

(A-35) 
(A-36) 
(A-37) 
(A-38) 
(A-39) 
(A-40) 
(A-41) 



(A-42) 
(A-43) 



(A-44) 
(A-45) 

first order of 



AV U3 = 



M di A i3 + A\M dj 



Ml 



AV 



Rij 



M di A\ + AjjMdj 
M 2 dl -M% 



for % ^ j, 



for i 7^ j. 



(A-46) 
(A-47) 



At this order M d is the physical down-type quark mass matrix, since there are no corrections 
to diagonal terms. 
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We derive relationship between Vq KM and Vckm- Since the W boson coupling is given by 

92 rrr tta+. 



C w ± = ^=[U L W+YD L + h.c] 

= ^=[u' l w;yvS km Vl(d)d' l + h.c] 



we obtain 



Vckm = ^ckm^l(-D)- 
Explicit form for each elements are given by 



^CKMll 


— Vckm 11, 


V CKM12 


= Vckm 12 j 


^CKMIS 


— Vckm 13 


,/ CKM21 


= Vckm 21 , 


> / CKM22 


= Vckm 22; 


K CKM23 


— Vckm 23 



1 + [E^+E~ h ^]tan(3 



1 + [^0) + ^(33) + £,(13) + ^(13) + £, (133)] ^ ^ 



1 + [^ (3) +£r (33) ]tan/3 



^(33) + £,(23) + ^(23) + £, (233)] ^ 



K CKM31 



V CKM32 



Vckm 



= Vckm 



32" 



l + [E/ 3) + ^ (33) ]tan/? 
1 + [^ (3) 
l + [E~W + E~ h ^]tan(3 



EL 31 ] tan (3 



K CKM33 



l + [^ 3 )-£f 2 )]tan/?' 

Vckm 33- 

Then the charged Higgs couplings can be expressed as 

C H ± = smpH~D Ri y dij U Lj + h.c. 

= smPH-D^VliD^yakV^U)^ + h.c. 
= smPH-D R 'yi(D) tk y dk V°i Mk3 U L ' 3 + h.c. 

Using Eqs.(A-47) and (A-51) - (A-58), we obtain 

M dl 



C H ± fa — taxiBH D Ri — 
H v M *l + [£fc«]tan/3 



(A-48) 
(A-49) 

(A-50) 
(A-51) 

(A-52) 

(A-53) 
(A-54) 

(A-55) 

(A-56) 

(A-57) 
(A-58) 



(A-59) 



+ h.c. for = (1, 1), (1, 2), (2, 1), (2, 2), 

(A-60) 
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C H± « ^ t^pH-D^ i - ^ gjL_ — 4m,^l; + h.c. for = (3,1), (3, 2), 

(A-61) 

r ^zr-n ' ^* i + [ jg J3 )+jg j33) ]tan/3 

H± * t; ^ Rl 1 + EM tan/5 1 + [E-^ + ^ + ^ + E£*> + EL™] tan/3 ^ 

+h.c. for (i,j) = (l,3),(2,3), (A-62) 

^ - V tan^-D R ; i - [E j t) ^ — 4m,^l; + h.c. for (i, j) = (3, 3). (A-63) 

Notice that b — > c and 6 — > u transition do not receive the higgsino loop contribution just as 
in the simplest case discussed in section 2. 
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